To appear in The Astrophysical Journal 

Preprint typeset using L^T^X style cmulatcapj v. 04/20/08 



ON 
O 
O 

> 

O 



< 
O 

6 
to 



BLAST OBSERVATIONS OF RESOLVED GALAXIES: TEMPERATURE PROFILES AND THE EFFECT OF 

AGN ON FIR TO SUBMILLIMETER EMISSION 

Donald V. Wiebe, 1,2 Peter A. R. Ade, 3 James J. Bock, 4,5 Edward L. Chapin, 1 Mark J. Devlin, 6 Simon Dicker, 6 
Matthew Griffin, 3 Joshua O. Gundersen, 7 Mark Halpern, 1 Peter C. Hargrave, 3 David H. Hughes, 8 
Jeff Klein, 6 Gaelen Marsden, 1 Peter G. Martin, 91 " Philip Mauskopf, 3 Calvin B. Netterfield, 2,10 Luca Olmi, 11 ' 12 
Enzo Pascale, 3 Guillaume Patanchon, 13 Marie Rex, 6 Douglas Scott, 1 Christopher Semisch, 6 Nicholas Thomas, 7 
Matthew D. P. Truch, Carole Tucker, 3 Gregory S. Tucker, 14 Marco P. Viero 10 

To appear in The Astrophysical Journal 

ABSTRACT 

Over the course of two nights, the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) 
made resolved maps of seven nearby (< 25Mpc) galaxies at 250, 350, and 500 p. During its June 

2005 flight from Sweden, BLAST observed a single nearby galaxy, NGC 4565. During the December 

2006 flight from Antarctica, BLAST observed the nearby galaxies NGC 1097, NGC 1291, NGC 1365, 
NGC 1512, NGC 1566, and NGC 1808. We fit physical dust models to a combination of BLAST 
observations and other available data for the galaxies observed by Spitzer. We fit a modified blackbody 
to the remaining galaxies to obtain total dust mass and mean dust temperature. For the four galaxies 
with Spitzer data, we also produce maps and radial profiles of dust column density and temperature. 
We measure the fraction of BLAST detected flux originating from the central cores of these galaxies 
and use this to calculate a "core fraction," an upper limit on the "AGN fraction" of these galaxies. 
We also find our resolved observations of these galaxies give a dust mass estimate 5-19 times larger 
than an unresolved observations would predict. Finally, we are able to use these data to derive a 
value for the dust mass absorption co-efficient of k = 0.29 ± 0.03 m 2 kg -1 at 250 pm. This study is 
an introduction to future higher-resolution and higher-sensitivity studies to be conducted by Herschel 
and SCUBA-2. 

Subject headings: balloons — galaxies: photometry — submillimeter — telescopes 
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1. INTRODUCTION 

Much of the bolometric luminosity produced by en- 
tire galaxies is emitted as thermal radiation from dust 
grains in the interstellar medium (ISM) at temperatures 
of several tens of Kelvin. The heat source is predomi- 
nantly optical and ultra-violet (UV) light from stars, as 
well as matter accretion onto an active galactic nucleus 
(AGN) in some objects. In either case, the spectrum of 
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this light tends to peak in the range ^60-200 um. Since 
the launch of the Infrared Astronomical Satellite (IRAS) 
in 1983, observations at these wavelengths have been es- 
sential in a number of fields, ranging from the studies of 
the earliest stages of star-formation in dusty molecular 
clouds in our own Galaxy, to the detection of distant star- 
forming galaxies that emit most of their light in this pre- 
viously unexploited band. Furthermore, excluding light 
from the Cosmic Microwave Background, it is known that 
approximately half of all of the light in the Universe is 
emitted at thes e wavelengths, the Cosmic Infrared B ack- 
gro und fCIB) (iPuget et alj |1996t iFixsen et all 1 19981 see 
also lDevlin et al.ll2009l) ~ 

Observations from the ground are impossible at wave- 
lengths in the range ~60-200 um, the far infrared (FIR) 
band, due primarily to the absorption and emission 
from atmospheric water vapor. Therefore, most of the 
progress at these wavelengths since IRAS has also been 
undertaken by satellites: the Infrared Space Observatory 
(ISO), the Spitzer Space Telescope, and the Akari Space 
Telescope. However, the lack of data at longer wave- 
lengths can bias the interpretation of these results to- 
ward the assumption of warmer dust temperatures, re- 
sulting in underestimates of luminosity and mass. At 
slightly longer wavelengths, ^350-1000 um, the submil- 
limeter band, several spectral windows are available 
from dry, high-altitude sites such as Mauna Kea or 
the Atacama Desert. In particular, studies at 850 um 



with SCUBA ( e.g. ISmail et alj|1 997: Hughes et alj|1998t 
iHolland et"aL1 Il999t lEales et alj 120001 : iDunne k Eales 



pOOttlScott et a l. 2002; Mort ier et al.ll2005l : IStevens et all 
\200d . iNutter fc Ward-Thompson! I2007h have been used 
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Fig. 1.— BLAST (bottom row) and MIPS (top row) observations of NGC 1097. The radius of the clipped BLAST map is 6'. The 
intensity scale is logarithmic. The FWHM size of the BLAST and MIPS beams are plotted for reference in the lower right corner. BLAST 
resolves the spiral arms and central bar of the galaxy. The 24 um image is shown for comparison only, and not used in this analysis. 



to detect very cool dust, both in our Galaxy, and in more 
distant galaxies. Typically in these surveys a dust tem- 
perature must be assumed in order to estimate redshifts 
and luminosities, because data are available only at a 
single wavelength. Therefore the combination of ground- 
based submillimeter observations with space-based FIR 
measurements has been required to assess the bolomet- 
ric luminosities and apparent dust temperatures of ev- 
erything from individual pre-stellar objects, up to entire 
galaxies at redshifts z > 1. However, the general lack of 
high signal-to-noise data in the crucial wavelength range 
~200-500um results in large uncertainties in these quan- 
tities. 

In this paper we present results from a study at 250, 
350, and 500 pm of several nearby galaxies using the 2-m 
Balloon-borne Large Aperture Submillimeter Telescope 
(BLAST). BLAST was built to fill the wavelength gap 
between ground-based submillimeter and space-based 
FIR instrumentation , using a prototyp e of the SPIRE 
camera for Herschel (jGriffm et alj[2007h . 

In June 2005, BLAST was launched from Esrange, 
Sweden, for it's first scientific flight, hereafter BLAST05. 
It acquired data for 100 hours before landing on Vic- 
toria Island, Canada. During this flight, BLAST ob- 
served the nearby galaxy NGC 4565. In Decem- 
ber 2006, BLAST conducted its second science flight, 
BLAST06, carrying out 250 hours of observations above 
Antarctica. Six nearby galaxies were observed dur- 
ing BLAST06: NGC 1097, NGC 1291, NGC 1365, 



NGC 1512, NGC 1566, and NGC 1808. 

All these galaxies are resolved by BLAST. The full 
widths at half-maximum (FWHM) of the BLAST06 
beams are 36, 42, and 60" at 250, 350, and 500 pm 
respectively. The BLAST05 fli ght suffered from degra- 
dation of the warm optics (jTruch et al.l 120081 ). The 
single galaxy in our sample from BLAST05 is never- 
theless resolved, despite the large, non-Gaussian point 
spread functions due to the degraded optics during this 
flight. We are therefore able to map the temperature 
and luminosity distributions, and hence dust column 
density, with greater signal-to-noise ratio (SNR) than 
any existing work. Our study can be cons idered a pre- 
curso r to The Herschel Reference Survey (|Boselli et al.l 
12008ft and the JCM T Nearby Galaxies Legacy Survey 
(jWilson et al.ll200"9T ) that will be undertaken by Herschel 
and SCUBA-2, respectively. A co mplete description o f 
the BLAST instrument is given in iPascale et aTT(|2008ft . 

The known correlation bet ween SEP shape and lumi- 
nosity in IRAS galaxies (e.g. iSoifer fc Neugebauerlll99lft 
means good measurements of these quantities for rest- 
frame galaxies are required to put proper constraints on 
high-redshift FIR source count models (e.g. Chapin et 
al., in prep.). BLAST samples the Rayleigh- Jeans side 
of thermal emission from local galaxies; when combined 
with FIR measurements from other instruments, the dust 
thermal emission peak is bracketed, allowing robust de- 
termination of these quantities. 

Resolved studies of galaxies are required to better 
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Fig. 2.— BLAST (bottom row) and MIPS (top row) observations of NGC 1291. The radius of the clipped map BLAST map is T. The 
intensity scale is linear. The size of the BLAST and MIPS beams are plotted for reference in the lower right corner. The outer 4' ring of 
spiral arm structure is detected in all three BLAST bands, as is the central core. 



model galactic structure and dynamics. The unprece- 
dented ability of BLAST to make high resolution maps 
on large angular scales in the submillimeter band pro- 
vides the best opportunity to-date to study the struc- 
ture of diffuse dust in nearby galaxies. To properly es- 
timate star formation rates from submillimeter fluxes of 
unresolved, high redshift galaxies, an understanding of 
the fraction of submillimeter flux originating from AGN- 
driven dust heating is required. Although BLAST is un- 
able to distinguish AGN-driven dust heating from nu- 
clear starburst driven heating, and this study involves 
nearby spirals, rather than high-redshift starbursts, we 
are still able to put limits on this AGN fraction by com- 
paring flux measured from the core of these resolved 
galaxies with their total integrated flux. Furthermore, 
because BLAST can make high-fidelity large-area maps, 
it is able to measure the complete emission from these 
spatially extended galaxies. 

2. DATA REDUCTION 

The common data redu ction pipelin e for B LAST 
has been presented in iTruch et al.l (|2008f ) and 
iPatanchon et al.l f2Q08i). Maps used in this analysis for 
NGC 4565 were ma de using the maximum -likelihood 
map maker SANEPIC (jPatanchon et all l2008| ) which was 
written explicitly for BLAST analysis. 

Maps of BLAST observations of the remaining g alax- 
ies were produced using ALMAGEST (|Wiebd I2008T) . an 
iterative, maximum-likelihood map maker also writ- 



ten for BLAST analysis, and algorithmically similar 
to o ther iterative map makers used in CMB analysis 
(e.g. lPrunet et a 1. 1 [2( )()()■. Because most of these galaxies 
were observed only at a single scan direction, and hence 
lack cross-linking, the maps produced are artificially con- 
strained to zero flux far from the observed galaxy to im- 
prove conve rgence properties of th e iterative algorithm 
(see §3.8 of IPatanchon et aH I2008D . In relatively high 
SNR cases, such as these galaxy maps, ALMAGE ST out- 
put i s consistent with maps made with SANEPIC (jWiebd 
HH), but may be produced faster, and allow arbitrary 
constraints to be applied simply. 

3. OBSERVATIONS 

3.1. BLAST Observations 

The seven nearby galaxies observed by BLAST were 
selected based on a number of competing criteria. First, 
the visibility requirements were restrictive, as they took 
into account BLAST'S limited elevation range (25°-60°), 
Sun and Moon avoidance criteria, and still had to make 
allowance for the changing latitude and longitude of the 
telescope arising from BLAST'S uncertain ground path 
while aloft. Second, the galaxies had to be sufficiently 
large to be resolved by BLAST'S optics (e.g. BLAST'S 
36" beam size at 250 pm corresponds to ~2.5 kpc at 
the average distance of 14 Mpc for this galaxy sample). 
Third, they had to be sufficiently luminous for BLAST to 
be able to map them with a high SNR in the relatively 
small fraction of the flights scheduled for these obser- 
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Fig. 3. — BLAST observations of NGC 1365. Pixels are 9" on a side. The radius of the clipped map is 6.'5. The intensity scale is 
logarithmic. 



vations. Finally, the sample needed to have measure- 
ments from the Infrared Astronomical Satellite (IRAS), 
the Spitzer Space Telescope, the Submillimeter Common- 
User Bolometer Array (SCUBA), or other experiments 
complementary to BLAST. 

Table [1] summarizes relevant physical parameters of the 
BLAST nearby galaxy sample* 5 The measured total flux 
density in each BLAST band is tabulated in Table [2] 

3.2. Other Multiwavelength Data 

As part of the Spitzer Infrared N earby Galaxies Sur- 
vey (SINGS: iKennicutt et al.l 120031 ). four of our sample 
(NGC 1097, NGC 1291, NGC 1512, and NGC 1566, 
hereafter BLAST/SINGS galaxies), have also been ob- 
served by the Multiband Imaging Pho tometer for Spitzer 
(MIPS) at 160 urn, 70 pm, and 24 urn (|Rieke et al.ll2004D 
and the Infrared Arra y Camera ( I RAC) at 8 pm, 6 pm, 
4.5 pm, and 3.6 pm jFazio et all 120041 ). both on the 
Spitzer Space Telescope. The remaining galaxies in our 
sample were not observed by Spitz er. Spatial l y inte - 
grated Spitzer fluxes are taken from IDale et all pOOl . 
and MIPS maps used in this analysis come from the fifth 
SINGS public data released 

IRAS data for all of our sample except NGC 1291 have 
been extracted from SCANPU 7 IRAS data for the fainter 
and more extend ed NGC 1291 are me asured from HIRES 
generated maps (|Aumann et all "l990). These IRAS data 
are used, in addition to Spitzer, to constrain the spatially 
integrated SEDs of these galaxies. A 20% uncorrelated 
uncertainty is applied to all IRAS data. IRAS 25 pm 
data are not used in this analysis when equivalent MIPS 
24 pm data are available. 

Two of our sample have been obser ved by SCUBA: 
NGC 1097 at 850 pm (IDale et all 12005ft. and N GC 1808 
at both 450 and 850 pm ([Stevens et al.l 120051 ). These 
measurements do not cover the entire region of emission 
detected by BLAST and we omit these measurements 

15 Instructions for obtaining BLAST maps used in this analysis 
may be found on http://blastexperiment.info/ 

16 obtained from http://data.spitzer.caltech.edu/popular/ 
sings/ 

17 available at http : //scanpi . ipac . caltech. edu: 9000/ 
appl i cat ions/Scanpi/ index .html 



from our fits for this reason. 

F or NGC 1808, we ca lculate aperture corrections for 
the lSteyens et al.l (2 005) reported flux densities. We ap- 
ply tdiel^eyen^^^Oapertures to the BLAST observa- 
tions and compare the 500 pm flux inside these apertures 
to the total 500 pm flux. From this, we calculate flux den- 
sities for NGC 1808 of s 450 = 13.12 Jy and s 850 = 1.76 Jy 
by applying the aperture corrections to the SCUBA mea- 
surements. These values are consistent with our (3 = 2 
modified blackbody model fit to BLAST plus 100 pm 
IRAS (SI4.3p. The aperture-corrected SCUBA measure- 
ments are plotted in Fig. [8l along with the model fit. 

We also exclude data from the Infrared Space Observa- 
tory (ISO) for NGC 1365 and NGC 4565, again because 
of aperture effects. 

3.3. Individual Galaxies 
A brief description of each observed galaxy follows. 

3.3.1. NGC 1097 

Features of NGC 1097, presented in Fig. [TJ include 
a bar 20 kpc in length that tapers off into two spiral 
arms (Bec k et al.l 120051) , a s well as several fringe struc- 
tures called jets or 'rays' (|Wolstencroft fc Zcaley 119751 : 
Eot71[]I1|. NGC 1097 has an active Seyfert 1/LINER 
nucleus, which is c onsidered to have a low luminosity 
(jMason et al.1 12007). and is surrounded by a 1.5 kpc di- 
ameter n uclear ring, contai ning 1.1 x 10 9 M^ of molec- 
ular gas (jHsieh et al.l 120081 : iGerin et al.l ll988l Thermal 
dust h eating in the nucle ar ring is dominated by a star- 
burst (M ason et al.ll2007t) with a st ellar production rate 
of SMgyr" 1 (jHummel et al.l 119871 ). BLAST made two 
observations of a ~0.25 deg 2 area centered on NGC 1097, 
with a total integration time of 52 minutes. BLAST 
clearly resolves the bar and spiral arms of the galaxy, 
but has insufficient resolution to resolve the 1.5 kpc star 
forming ring. 

3.3.2. NGC 1291 

NGC 1291, presented in Fig. H is t he brightest Sa 
galax y in the Shapley-Ames Catalogue (Brcgman et al. 
1995). NGC 1291 has an orientation that is almost 'face 
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Fig. 4.— BLAST (bottom row) and MIPS (top row) observations of NGC 1512. The radius of the clipped BLAST maps is 6'. The 
intensity scale is linear. The size of the BLAST and MIPS beams are plotted for reference in the lower right corner. BLAST detects the 
companion galaxy, NGC 1510 (a = 4 h 03 m 33 s , 6 = -43°24'; lower right) at 250um and 350 urn. In these two bands, BLAST also can 
distinguish between the core and disk of NGC 1512. 



on' and has a primary bar and second ary bar that have 
a 30° misalignment with each other ([Perez fc Freeman! 
I2006D . The total HI mass is 2.4 x 10 9 M Q , and there 
is little evidence fo r star formation in the central bulge 
(Hogg et al.ll200lT) . BLAST made four observations of a 
~0.4deg 2 area centered on NGC 1291, with a total in- 
tegration time of 92 minutes. Although diffuse, BLAST 
detects both the central core of the galaxy, as well as the 
ring-like spiral structure 4' from the core. 

3.3.3. NGC 1365 

NGC 1365, presented in Fig. [3l is is usually classified 
as a Seyfert 1.5, bu t some authors list it as Seyfert 1 
or Se yfert 2 (e.g., [Galliano et al] 120051 : iVeron et all 
[19801: lAlloin et all 119811 ). The AGN has low luminos- 
ity, but there is evidence from optical studies to sug- 
gest that there is star form ation surrounding the nucleus 
(jKomossa fc Schulzl[l998f ). BLAST made three observa- 
tions of a ~0.25deg 2 area centered on NGC 1365, with 
a total integration time of 30 minutes. The brightest 
of our sample, BLAST'S observations prominently show 
both the bar and spiral arm structure of the galaxy. 



3.3.4. NGC 1512 



has a nucleus sur- 



NGC 1512, presented in Fig. 
rounded by a circumnuclear ring (Ma oz et al.|[2000T ). Be- 
yond the ring, its spiral arms are complicated by its 



satellite galaxy, NGC 1510 (|Kinmanl [1978?) . BLAST 
made three observations of a ~0.45 deg 2 area centered 
on NGC 1512, with a total integration time of 28 min- 
utes. BLAST detects the core and tightly wound central 
spiral of NGC 1512. BLAST also detects the companion 
galaxy, NGC 1510. 

3.3.5. NGC 1566 

NGC 1566, presente d in Fig. El is the seco nd brightest 
known Seyfert galaxy (jde Vaucouleursl[l973h . Extending 
from the nucleus are two prominent spiral arms which 
continue outward to form the gal axy's outer ring that 
is at an inclin ation of about 30° (jde Vaucouleurslll975l : 
lBottemafL992T ). The arms cont ain a significant amount of 
star formation (|Bottemalll992l ). BLAST made two obser- 
vations of a ^0.7 deg 2 area centered on NGC 1566, with 
a total integration time of 50 minutes. BLAST resolves 
the spiral structure of the galaxy, and the observations 
are dominated by submillimeter flux originating in the 
core of the galaxy. 

3.3.6. NGC 1808 

NGC 1808, presented in F ig. El is a Seyfert 2 galaxy 
(| Jimenez- Bailon et al.l 12005) . Subsequent measurement 
seems to indicate that a large part of the radiation ema- 
nating from the central 1 kpc of the galaxy is from an 
active gala ctic nucleus as well as a high rate of star 
formation (|Forster Schrcibe r et al.l [20031 : Maiolin o et al.l 
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Fig. 5.— BLAST (bottom row) and MIPS (top row) observations of NGC 1566. The radius of the clipped BLAST maps is 5'. The 
intensity scale is linear. The size of the BLAST and MIPS beams are plotted for reference in the lower right corner. BLAST detects the 
spiral arms, although the map is dominated by flux from the core. 

at 250 pm, 350 pm, and 500 pm, respectively. The 250 pm 
fit, which has the highest signal-to-noise ratio, agrees well 
with scale lengths measured at other wavelengths (e.g. , 
van der Kruit & Seari3l981tlEngargiola & Harpeifll992t 
Rice et al 119961 ). In this fit we have ignored the central 3' 
region of the galaxy. This region corresponds to an area 
of reduced 250 pm flux density, po ssibly corresponding t o 
the molecular ring seen in 12 CO (INeininger et al.lfl999h 

4. SED FITTING 

4.1. Physical Dust Models 

For the four BLAST/SINGS galaxi es, we fit 
the si licate-graphite-PAH dust models of iDraine fc Lil 
(2007), hereafter DL07, to BLAST, MIPS, IRAS, and 
IRAC data. These models have already been applied, 
without t he benefit of BLAST observations, to the same 
galaxies in lDraine et al.l (|2007l ). hereafter Dr07. We inde- 
pendently recreate this analysis with the inclusion of the 
BLAST data. These models provide the dust emissivity 
per hydrogen atom, ju(qpAMiU mia ,U max ), dependent on 
three model parameters: 

1- Qpah, the fraction of total dust mass in poly- 
cyclic aromatic hydrocarbons (PAHs) containing 
less than 10 3 carbon atoms; 

2. [/ m i n , a starlight intensity factor characterizing the 
radiation field heating the diffuse ISM; 



l200l . Using ROSAT, Uunkes et all (fl995h estimate 
a star formation rate of between 5 and 13MQyr _1 , 
and a supernova event rate of 0.09 yr _1 . A companion 
galaxy, not observed by BLAST, NGC 1792 is located 
40' from NGC 1808 (130 kpc at a distance of 10.9 Mpc) 
and may be respon sible for an accelerated star formation 
rate in NGC 1808 (jJimenez-Bailon et al.|[2005l ). BLAST 
made four observations of a ~1 deg 2 area centered on 
NGC 1808, with a total integration time of 179 min- 
utes. BLAST detects both the core and disk of the 
galaxy. BLAST also detects excess flux in this map near 
the location of the interacting galaxies ESO 305-IG 010 
(a = 5 h 08 m 27 s , 5 = -37°39f5, z = 0.0524). 

3.3.7. NGC 4565 

NGC 4565, presented in Fig. is an edge-on barred- 
spiral galaxy in the constellation Coma Berenices. It has 
an apparent size of about 12' x 2'.5. Although degraded 
by the large point-sprea d functions of BLAST05 (186- 
189", iTruch et al.l[2008h . NGC 4565 is still resolved by 
BLAST. BLAST made one observation of a ~0.4 deg 2 
area centered on NGC 4565, with a total integration 
time of 49 minutes. A beam-d econvolved image, us ing 
the method described in §2.4 of iChapin et al.l (|2008h . is 
shown in the bottom row of Fig. [7J 

We fit an exponential profile to the beam-deconvolved 
maps of the form exp(|a;|/s), where a; is a spatial variable 
aligned along the major axis of the galactic disk and s is 
the scale length, and find that s = 118", 156", and 142" 
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Fig. 6. — BLAST observations of NGC 1808 (upper right). Pixels are 9" on a side. BLAST also detects excess flux density near 
the position of the interacting galaxies ESO 305-IG 010 (a = 5 h 08 m 27 s , 6 = -37°39.'5, z = 0.0524; lower left). The intensity scale is 
logarithmic. 



3- t/ max , a starlight intensity factor characterizing the 
radiation field in more intense star forming regions, 
including photodissociative regions (PDRs). 

Both U m i n and U max are relative to the specific energ y 
density of starlight measured by iMathis et al.l (|1983f) . 
Dr07 consider a fourth model parameter, a, which is 
the exponent characterizing the power-law decrease in 
starlight intensity: dM d /dU cx U~ a . They find the 
model fits to be insensitive to its value, and the mod- 
els provided have a — 2. 

Our fitting procedure follows that of Dr07. We re- 
strict ourselves to the seven M ilky Way dust model sets 18 
(jWeingartner k. Drainell2001l updated by DL07). These 
seven model sets have values of (?pah between 0.47% and 
4.6%. We take linear combinations of these to produce 43 
different model sets for a finer-grained sampling of (/pah- 
Each set contains values of j v for twenty-two values of 
t^min between 0.1 and 25, and for five "PDR models" 
with values of J7 max (10 2 , 10 3 , 10 4 , 10 5 , 10 6 ), plus one 
"diffuse ISM model" with U max = U m i n . 

These dust models are combined with a stellar black- 
body to produce a galactic model: 



-^1/, model — Q+B U (T+ 



M ri 



Mh 
M d 



X [(! - l)ju{QpAM, U min , U min ) 
+7>(?PAH, U mm , E/max)] , 



(1) 



where fi* is the solid angle subtended by stars, B y {T) 
is the Planck function, the nominal stellar temperature 
T± = 5000 K, Md is the total dust mass, tor the mass of 
a hydrogen atom, D the distance to the galaxy, and 7 a 
mixing fraction between the diffuse ISM (L7 max = £/min) 
and PDR ((7 max > U m i n ) dust models. The hydrogen-to- 
dust mass ratio, Mn/Md, is set by the model choice, and 
varies from 96 to 100. 
Like Dr07, we find that the best fits for the four galax- 



ies are insensitive to the choice of U n 



and so set 



18 obtained from http://www.astro.prlnceton.edu/~dralne/ 
dust/ dust . html 



CAnax = 10 6 . We obtain the best-fit values of fi*, Md, 
7, g PA H, and U m - m through \ 2 minimization (see § 14. 5|) . 

The best-fit models for the galaxies to the extended 
Spitzer- IRAS '-BLAST data set are plotted in Fig. [H 
This figure also plots the spatially integrated Spitzer 
and IRAS, BLAST, and SCUBA measurements for these 
galaxies. 

The model parameters are tabulated in Table El We 
also tabulate the dust-weighted mean starlight intensity 
scale factor 



(U) 



(1 - 7)C/min + 7 



ln([/max/[/ mi n) 



u min (-'max 



(2) 



and the dust luminosity calculated from the model fit. 

4.2. Modified Blackbody SEDs 

In addition to the DL07 models, we also fit BLAST 
plus available MIPS 70 and 160 pm and IRAS 100 pm 
data to a single-component modified blackbody SED: 



odd 



M d K 







B V {T), 



(3) 



where the dust emissivity index, (3, is fixed to 2, and vq 
to 1.2 THz = c/250um. An uncertainty of ±0.3 Mpc is 
assumed for all galaxy distances. This simplified SED 
model is needed in the absence of Spitzer data, where 
the DL07 models are insufficiently constrained. We also 
use this model in our subgalactic SED fits fS I4.4p due to 
its computational simplicity. 

The dust mass absorption co-efficient at vo, k, is not 
well determined, and depen ds on the object investigated 
(see iNetterfield et al.l 120091 for a discussion) . For each 
BLAST/SINGS galaxy, we calculate a value for k by 
comparing the DL07 fit to the fit found for the spa- 
tially integrated modified blackbody SED (§[4~3]). The 
values of k for each galaxy are summarized in Table [3l 
These values are used in the calculation of the column 
density maps presented in Fig. [TO] For the galaxies 
without Spitzer observations, we use the mean value, 
K 250 = 0.29 ± 0.03 m 2 kg -1 , to calculate dust mass 
from the modified blackbody fits. Assuming a n oc v 2 
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Fig. 7. — Top row: BLAST observ ations of NGC 4565. Bottom row: Results of deconvolving the raw maps following the procedure 
outlined in §2.4 of Chapin ct al. (2008). This deconvolution process results in 3—5 times better resolution, at a cost of reduction in signal- 
to-noise of a factor of ~2. Pixels are 20" on a side . Note the linear intensity scale. Effective beam sizes, based on the full widths at 
half-power (FWHP) found in Table 1 of Chapin ct al. (2008]), are pl otted in the lower r ight corners for both the raw and deconvolved maps. 
The full BLAST05 point spread functions are shown in Figure 1 of l'lruch et al.l 1120081 ). 

For the four BLAST/SINGS galaxies, we also fit Equa- 
tion [3] to BLAST, plus 160, 100, and 70 pm data from 
MIPS or IRAS, to compare with the results of the DL07 
models. These fits are plotted in Fig.[8jas dashed curves. 
We note reasonably good agreement between the two 
models for all galaxies near the thermal peak. On the 
Wien side of the peak, as expected, the single-component 
modified blackbody quickly falls off in comparison to the 
DL07 model fit. On the Raleigh- Jeans side, due to the 
adoption of an emissivity index (3 = 2, the modified 
blackbody also falls off faster than the DL07 model fit. 
For NGC 1097 and NGC 1291, the modified blackbody 
fit underpredicts the 500 pm flux density by more than 

1-CT. 

Dust temperatures from these fits for all galaxies in the 
sample are listed in Table [21 as are dust masses for the 
three galaxies without Spitzer observations. 



evolution, this corresponds to k S50 = 0.025 m 2 kg -1 , 
sm aller than the fcssn = 0.07 m 2 kg -1 value calculated 
by I James et al.l (I2002D. It is a lso l ower than the 250 pm 
values of IDraine fc Led (fl98l and lCasevi (|l99ll ). 

For each fit, a brute- force grid search is performed to 
simultaneously fit and T to the data. The temper- 
ature, T, is sampled at 1000 uniformly-spaced points 
between 5K and 30 K. This range was determined by 
first letting temperatures vary over a much wider range 
(1K-100K) and then discarding those portions of the 
parameter space discovered to be irrelevant to the fit- 
ting procedure. Because the dust mass, M<j, may vary 
by several orders of magnitude, the sampling limits for 
the dust mass are automatically adjusted as necessary 
for the fit, and sampled uniformly at 1000 points. As 
with the DL07 models, goodness-of-fit is assessed using 
a x 2 test (§ !4.5p . The marginalized likelihood is used to 
determine medians and 68% confidence intervals for the 
two fit parameters. 

4.3. Spatially Integrated SED Fits 

For the three galaxies without Spitzer data, we fit 
Equation[3]to BLAST plus IRAS 100 pm. The IRAS and 
BLAST measurements used to constrain the fits, and the 
results from these fits are presented in Fig. [5J which also 
plots additional data from IRAS, ISO and SCUBA, not 
used in this analysis. 



4.4. Subgalactic SED Fits 

Because the MIPS beams at 70 pm and 160 pm are of 
comparable size to the BLAST beams, we also investigate 
the dust properties on sub-galactic scales for the four 
BLAST/SINGS galaxies. To do this, the BLAST maps 
produced for the four galaxies were made on the same 9" 
pixel grid as the MIPS 160 pm maps. The 70 pm MIPS 
images, which have a compatibly aligned 4'.'5 pixel grid, 
arc rcbinned to 9" to put them on the same grid. All five 
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Fig. 8. — Spatially integrated SEDs for the seven galaxies in our sample. Data from BLAST, MIPS, IRAC, and IRAS used in this analysis 
are shown as bl ack squares. Othe r data, which are n ot used in this analy sis, from IRAS for NGC 1365, NGC 1808, an d NGC 4565, SCUBA 
for NGC 1097 IDale et al.ll2005h , and NGC 1808 JStevens et al.l 120051) . and ISO for NGC 1365, and NGC 4565 ( Spinogli o~etaII |2002| ; 
IStickel et al.1120041) are the lighter diamonds. The aperture corrections ca lculated in jjl3T2lhav e been applie d to the SCUBA measurements 
of NGC 1808. The solid curves are the best-fit physical dust models of IDraine fc Lil 1 120071) (DL07; see §|4T). The dotted curve in the 
panel for NGC 1291 is the b est-fi t DL07 model subject to the constraint t/ m j n > 0.7. The dashed curves are the best-fit single-component 
modified blackbody model (§ 14.31 1, Parameter values for these fits are given in Tables [2] and [3] 
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maps are then convolved (smoothed) to the resolution 
of the 500 pm data (60") and our modified blackbody 
template (Equation [3]) is fit to the data at each pixel in 
the map. 

The central twelve pixels of NGC 1097 have 70 pm 
surface brightness in excess of 200MJysr _1 . For these 
twelve pixels, the 70 pm data are excluded, due to pos- 
sible non-linearities} 9 and the model is fit to four photo- 
metric points only. 

Resultant temperature and column density maps for 
the four galaxies are presented in Figs. [9] and [TU1 
NGC 1510, the companion galaxy to NGC 1512, ap- 
pears quite prominently in these maps as seemingly hot- 
ter (>25 K) than NGC 1512. This is not believed to be 
a real effect, and is likely due to poor performance by 
the fitting routine on the data for NGC 1510, for which 
= 2 is probably not an appropriate choice. 

Radially- averaged temperature and column density 
profiles are plotted in Fig. [TTJ We note a general de- 
crease in both column density and dust temperature with 
increasing radius. Of the four galaxies, we notice two 
subgroups of two galaxies each. The brighter NGC 1097 
and NGC 1566, for which BLAST detects high intensity 
emission continuously from the core to disk of the galaxy, 
have relatively flat temperature profiles. NGC 1291 and 
NGC 1512, which BLAST detects with a more isolated 
core, fall off more steeply. We also potentially detect the 
1.5 kpc star-forming ring in NGC 1097, which causes the 
kink in the temperature profile. 

4.5. x 2 Calculation 

We use x 2 minimization to assess goodness of fit 
for both the DL07 models and the modified blackbody 
model. The BLAST calibration uncertainti es are highly 
corre l ated between observational bands (jTruch et alJ 
2008, 2009). Given observational data, Sd, and band- 
convolved model predictions, s, we compute 



X 2 = (s - s d ) 



s<D T C~ 



(s - Sd) , 



(4) 



where C _1 is the inverse data covariance matrix, and (-) T 
denotes the transpose. Given the calibration uncertainty 
correlation matrix, p, the data covariance matrix is: 

dj = [(*d)< + «)*] + PH > (5) 

where Cd is the measurement uncertainty, <r c the calibra- 
tion uncertainty, and <r m an uncertainty associated with 
the model. For the DL07 models, in order to compare 
with their results, we follow Dr07 and adopt o~ m = O.lsd- 
For the modified blackbody model, we use er m = 0. We 
assume that IRAS and Spitzer uncertainties are uncor- 
rected between bands, and set pij = 5{j for non-BLAST 
measurements . 

For the BLAST measurements, we set p to the value of 
the Pearson c orrelation m atrices computed as outlined in 
ITruch et all (|2008l I2009D . For both the calibration un- 
certainty (<7 C ) and the calibration uncertainty correlation 
matrix (p), we add to the values given in those papers 
an additional source of uncertainty associated with our 
ability to measure the BLAST band passes. We esti- 
mate the error in the bandpass measurements to be at 

19 http: //data. spitzer . caltech . edu/popular /sings/ 

20070410_enhanced_vl/Documents/sings_f if th_delivery_v2 .pdf 



most 5%. We additionally assume that this uncertainty 
is uncorrelated between bands. 

The magnitude of the effect of bandpass measurement 
error on the BLAST calibration uncertainties depends on 
the spectrum of the object observed. A source with the 
same spectrum as our calibr ators (e.g., f o r B LAST06, 
VY Canis Majoris, for which ITruch et alJ (|2009h find a 
best-fit single-component modified blackbody with T = 
346 ± 19 K and = 0.55 ± 0.05) is unaffected by band- 
pass errors, but for a source with a significantly different 
spectrum, the effect can be non-negligible, especially in 
the high SNR case. This is the case in this study, and we 
have adjusted our calibration uncertainties accordingly. 

Because the underlying spectrum of the galaxies is un- 
known, we use the modified blackbody model as a proxy. 
The calculated calibration uncertainties affect the best 
fit of the model, so this is an iterative process: we first 
fit the model using a guess for the calibration uncertain- 
ties, and then update the uncertainties based on the best 
fit. This process is iterated until convergence is achieved. 
In general only one iteration after the initial guess is re- 
quired for convergence. 

For this procedure, we have calculated calibration un- 
certainties and uncertainty correlations for various = 2 
modified blackbody models with temperatures in the 
range 15-25 K. The resulting calibration uncertainties 
vary most with choice of model at 250 pm (12-14% for 
BLAST06), with smaller variations at 350 pm (10.0- 
10.5%) and 500 pm (8.00-8.25%). The uncertainty cor- 
relations calculated show similar behavior, varying be- 
tween 0.80 and 0.97 for the 250-350 pm correlation but 
only 0.91-0.92 for the 350-500 pm correlation. 

When reporting goodness of fit, we report x 2 per de- 
gree of freedom: 



1 



Xr 



N h -N r 



-x 



(6) 



where iVb is the number of observational bands consid- 
ered in the fit, and _/V m the number of model parame- 
ters (five for the DL07 models and two for the modified 
blackbody model). The number of bands ranges from 
four (3 BLAST + 100 pm IRAS) for modified blackbody 
fits in the absence of Spitzer data, to thirteen (4 IRAC, 
3 IRAS, 3 MIPS, 3 BLAST) when fitting the DL07 mod- 
els. 

5. DISCUSSION 

5.1. Dust Models 

Of the four galaxies observed by both BLAST and 
Spitzer, the DL07 fits we find for all but NGC 1291 are 
consistent with the best-fit models determined without 
BLAST data in Dr07. For NGC 1291 we find a dust 
mass approximately two times larger than Dr07. This 
discrepancy is due to the restriction U m i n > 0.7 imposed 
by Dr07 in the absence of SCUBA data. We fore-go 
this restriction for our fit of NGC 1291, and find a best- 
fit £/min = 0.2. Constraining U m i n > 0.7 results in a 
fit consistent with that of Dr07 (g P4H =2.8%, U m i n =Q.7, 
log(M d /M )=7.36, 7 = 1.00%), but a poor fit to BLAST 
data (x 2 = 2.59). This fit is plotted in Fig. [S] for com- 
parison. 

The dust temperatures found in this study agree 
well with the 15-25 K temperatures found for the 
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Fig. 9. — Dust temperature maps based on BLAST and MIPS observations for the four BLAST/SINGS galaxies. Only those pixels whose 
250 um flux density is inconsistent with zero at the 1-cr level are plotted. Contours from the 24 pm MIPS maps, smoothed with a 18"xl8" 
boxcar window function, are over-plotted for comparison. For NGC 1097 and NGC 1512, the O.SMJysr" 1 24-um contour is plotted. For 
NGC 1291, the O.lMJysr contour is used, and for NGC 1566, the 5MJysr — 1 contour is plotted. Mean temperature errors in the high 
SNR portions of the maps are ±0.8 K, ±2.5 K, ±1.1 K and ±0.8 K for NGC 1097, NGC 1291, NGC 1512, and NGC 1566, respectively. 
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5.2. Fraction of Core to Total Emission 

BLAST'S unique ability to spatially resolve the sub- 
millimeter flux of these nearby galaxies allows us to in- 
vestigate the fraction of submillimeter flux which arises 
from the extended disks compared to their central core 
and circumnuclear regions. Although BLAST'S resolu- 
tion is insufficient to resolve the central AGN of these 
galaxies (where present), we can still put limits on the 
contribution made to the submillimeter flux by the nu- 
clear accretion and star formation activity towards their 
circumnuclear regions (typically within a few kpc of the 
nucleus). 

For each galaxy, except the edge-on NGC 4565, we find 
the 3 db radial flux contour (that is, the radius at which 
the radially averaged flux has dropped to half of the peak 
flux) at 250 pm and use this radius to define the core of 
the galaxy. Fluxes in all three BLAST bands from this 
core region are presented in Table |U Also in this table 
is the "core fraction," the fraction of the total BLAST 
measured flux which originates in the core of the galaxy. 

Most of the sample have core fractions in the range 
7.5% to 15%. NGC 1566 is a clear outlier, with roughly 



one third of the submillimeter flux in the core of the 
galaxy. The second brightest known Seyfert galaxy 
(|de Vaucouleu7i[l973lL NGC 1566, has a core whose sub- 
millimeter flux originates largely from the AGN. The 
galaxy with the smallest core fraction is NGC 1291, a 
galaxy with no known active nucleus. Other than in 
NGC 1566, the bulk of the submillimeter emission in the 
core of the galaxy could be explained by nuclear star 
formation. 

All galaxies show a decrease in core fraction with in- 
creasing wavelength, a result of the core dust tempera- 
ture being generally warmer than the dust in the sur- 
rounding disk (see Fig. [TTj) . Other than NGC 1566, 
most galaxies have weak or no AGN, although BLAST 
is unable to separate AGN derived flux from star forma- 
tion in the central regions of the galaxies. In the case 
of NGC 1097, t he unresolved 18" circumnuclear star- 
burst ring (e.g.. IHummel et alj 119871: IBarth et al.lll995l : 
lOuillen et all 119951; iKotilainen et all (20001 ) . likely pro- 
vides a sizable contribution to the total core flux. 

BLAST'S ability to measure the compact to extended 
flux ratio of these galaxies is a result of BLAST'S high 
sensitivity and resolution. This indicates that cur- 
rent and future, higher-sensitivity and higher-resolution 
FIR and submillimeter experiments, such as Herschel, 
SCUBA-2, SPICA, JWST, ALMA, LMT, and GBT, 
which will be able to spatially separate non-thermal 
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Fig. 10. — Dust column density maps based on BLAST and MIPS observations for the four BLAST/SINGS galaxies. Only those 
pixels whose 250 um flux density is inconsistent with zero at the 1-cr level are plotted. Contours plotted are the same as in Fig. [9] Mean 
dust column density errors in the high SNR portions of the maps are ±2.1x10 Mq sr - 1 , ±1.6x10 13 Mq sr — 1 , ±0.9xl0 13 Mq sr — x , and 
±2.1 xlO 13 Mq sr" 1 for NGC 1097, NGC 1291, NGC 1512, and NGC 1566, respectively. 



1 Q sr 

components from extended star formation emission, are 
important for understanding the relationship between 
FIR-submillimetcr luminosity and star formation rates 
in both local and high-redshift galaxies. 

5.3. Comparison of Resolved and Unresolved Mass and 
Luminosity Estimates 

Our resolved measurements of the mass distribution in 
these galaxies allow us to investigate how mass estimates 
may be affected in more distant and smaller galaxies that 
are unresolved. For the four galaxies for which we have 
dust column density maps (Fig, mi)), we apply the same 
aperture used to calculate our total flux measurements 
to the column density maps. The masses contained in 
these apertures, along with the aperture diameters, are 
listed in Table El 

In all cases, the resolved mass estimate is significantly 
higher than the mass estimate obtained by the DL07 fit 
to the unresolved total flux from the galaxy. The ratio 
of resolved to unresolved mass estimates ranges from ~5, 
for NGC 1566, to -19, for NGC 1291. Furthermore, the 
mass ratio follows the behavior of the radial tempera- 
ture profile (Fig. Hip, with a larger temperature gradient 
corresponding to a larger mass ratio. This effect is inde- 
pendent of the choice of k, since both our resolved and 
unresolved mass estimates use the same k to derive dust 
mass from the SED fit. 

As a check, we performed a similar analysis to measure 



the difference between resolved and unresolved luminos- 
ity estimates. Using the same aperture, we summed the 
60-1000 pm luminosity in the resolved maps, and com- 
pared this to the 60-1000 pm found from the modified 
blackbody fit to the galaxy as a whole. As expected, 
these values agree with one another to within measure- 
ment uncertainty. 

Clearly, while luminosity, and hence star-formation 
rate, is unaffected by resolution, accurate dust mass mea- 
surements rely on resolved observations. Since lumi- 
nosity, and therefore observed flux density, is a strong 
function of temperature, there is a potential for unre- 
solved observations of galaxies to hide cool dust. For 
the four galaxies we investigate, a disproportionate frac- 
tion of their total SEDs is produced by their relatively 
warmer, less massive cores. A single temperature mod- 
ified blackbody fit to such data tends to infer a warmer 
temperature and shallower (i.e. smaller) value of 0, and 
hence lower, incorrect, dust mass. Our data are consis- 
tent with the model of iDunne fc Ealej (|200lD . who use 
a two-component (3 = 2 modified blackbody SED to sep- 
arate the cool (15-25 K) dust component from warmer 
(T > 30 K) emission in the SCUBA Local Universe 
Galaxy Survey (SLUGS). 

The upcoming higher-resolution experiments should be 
able to better estimate the quantity of dust in distant 
galaxies. 
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Fig. 11. — Mean radial dust temperature (top) and dust column 
density (bottom) for the four BLAST/SINGS galaxies, with 1-cr un- 
certainty bands. Resolution elements at 500 pm are approximately 
4.9, 5.2, 3.0, and 2.8 kpc for NGC 1097, NGC 1291, NGC 1512, 
and NGC 1566, respectively. 

6. SUMMARY 

BLAST made resolved observations of seven nearby 
galaxies over the course of two flights. Four of these 
galaxies, NGC 1097, NGC 1291, NGC 1512, and 
NGC 1566, have complementary Spitzer observations as 
part of the SIN GS survey. For the se four galaxies, we 
fit the models of lDraine fell (|2007fl to BLAST, Spitzer, 
and IRAS observations. Best-fit parameters are tabu- 
lated in Tabled 

The best-fit models for three of these four galax- 
ies are consistent wit h the best-fit models calculated 
in iDraine et ail (|2007f ) without BLAST data. For the 
fourth galaxy, NGC 1291, we find a dust mass roughly 
two times larger t han the dust mass determined by 
Draine et al. (2007). We also calculate a value for the 
dust mass absorption coefficient at 250 pm n = 0.29 ± 
0.03 m 2 kg -1 by comparing the Draine & Li models with 
a modified blackbody model. 

For these four galaxies, we also produce maps of dust 
column density and mean dust temperature based on 
BLAST and MIPS observations fit to a single compo- 
nent modified blackbody model (Figs. 191 and [TO")), as well 
as radial profiles of the same quantities (Fig. [TTj) . For the 
remaining three galaxies observed by BLAST but not ob- 
served by Spitzer, NGC 1365, NGC 1808, and NGC 4565, 



we calculate spatially integrated dust temperatures and 
dust masses by fitting a modified blackbody model to 
BLAST and IRAS data. We find mean dust tempera- 
tures for our galaxy sample to be in the range 16-23 K. 

We calculate the fraction of the submillimeter detected 
flux originating in the core of the galaxy as a fraction of 
the total submillimeter emission from the galaxy. Al- 
though BLAST is unable to resolve the nucleus of these 
galaxies, the ratio of compact to extended flux measured 
by BLAST puts an important upper limit on the fraction 
of submillimeter radiation driven by a central AGN. 

Only one of our sample, NGC 1365, the second bright- 
est known Seyfert galaxy, has a significant fraction of its 
flux (>30%) emanating from the core of the galaxy. In 
the remainder the core fraction is small (<15%) and may 
be accountable primarily to nuclear star formation. In all 
cases, the core fraction increases with decreasing wave- 
length, indicating that dust in the core of these galaxies 
is warmer than the dust in the surrounding disk, agreeing 
with the radial temperature profile of Fig. [TTJ 

We compare the total dust mass in our column density 
maps to the dust mass estimate based on the observa- 
tion of the galaxies as a whole (i.e., an unresolved obser- 
vation). We find 5-19 times more dust mass in the re- 
solved maps than the unresolved estimate, a result of the 
unresolved measurements hiding cool dust. We perform 
a similar analysis with luminosity and find no discrep- 
ancy between the resolved and unresolved 60-1000 pm 
luminosity estimates. 

Finally, for NGC 4565, we calculate a scale length of 
118", 156", 142" at 250pm, 350pm, and 500pm, respec- 
tively, in good agreement with measurements at other 
wavelengths. 
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BLAST Observations of Resolved Galaxies 



TABLE 1 

Physical parameters of BLAST nearby galaxies 



Galaxy 


R.A. 
(J2000) 


Dec. 
(J2000) 


Dist. 
[Mpc] 


TypeW 


Nucleus 


NGC 1097 


2 h 46 m 19 s 


-30°16'30" 


16.9( b > 


SB(s)b 


Syl 


NGC 1291 


3 h 17 m 18?6 


-41°06'29" 


9.7< b > 


(R)SB(s)0/a 




NGC 1365 


3 h 33 m 36?4 


-36°08'25" 


18.4< c ) 


SB(s)b 


Sy2 


NGC 1512 


4 h 03 m 54?3 


-43°20'56" 


10.4< b ) 


SB(r)a 




NGC 1566 


4 h 20 m 00?4 


-54°56'16" 


18.0< b ) 


SAB(s)bc 


Syl 


NGC 1808 


5 h 07 m 42 s 3 


-37°30'47" 


10.9< d ) 


(R)SAB(s)a 


Sy2 


NGC 4565 


12 h 36 m 20?8 


+25°59'16" 


12.6^ 


SA(s)b 


Sy2 



a RC3 morphological type from dc Vaucoulcurs ct al. (1991) 
b IKennicutt et al.l O(503T) 
c Komossa & Schulz (1998) 
d I Jimenez-Bailon et aLl I (|2005I 1 
c Lawrence ct al. (1999) 



TABLE 2 

Observed flux densities and derived quantities for BLAST nearby galaxies 



Galaxy 


250 pm 




350 pm 


500 pm 






[Jy] 




[Jy] 


[Jy] 


[K] 




NGC 1097 


89.4±0.2± 


7.7 


40.5±0.1±3.0 


16.2±0.1±1.4 


21.3±0.3 


... (a) 


NGC 1291 


26.1±0.2± 


2.6 


12.7±0.1±1.0 


6.4±0.1±0.5 


16.6±0.4 


... (a) 


NGC 1365 


145.8±0.3±12.9 


62.3±0.2±4.6 


24.7±0.1±2.1 


19.9±0.4 


a 7o +0.08 

s -'^-o.io 


NGC 1512 


14.7±0.1± 


1.3 


6.9±0.1±0.5 


3.4±0.1±0.3 


20.3±0.3 


... (a) 


NGC 1566 


51.5±0.2± 


1.6 


22.2±0.1±1.7 


9.1±0.1±0.8 


20.1±0.4 


... (a) 


NGC 1808 


84.0±0.1± 


7.1 


33.1±0.1±2.5 


11.8±0.1±1.0 


22.8±0.5 


7.851°;?? 


NGC 4565 


37.2±0.5± 


4.5 


21.0±0.3±2.1 


9.8±0.3±0.8 


16 1+ ' 4 
iD - i -0.5 


8.2llg;S 



Note. — Flux uncertainties reported are first the measurement error (<r c i), and 
then the calibration error (u c ). Temperatures and dust masses are extracted from the 
modified blackbody fits presented in Fig. [8] 



a This modified blackbody fit is used to calculate k from the DL07 model fit, so no 
independent mass is determined. 



TABLE 3 

DL07 model parameters for BLAST nearby galaxies 



Galaxy lc 






<?PAH 

% 


(U) 


^min 


7 
% 


xl 


K (a) 

m 2 kg -1 


NGC 1097 


8.37 


10.80 


3.0 


2.00 


1.5 


2.80 


0.58 


0.299 


NGC 1291 


7.62 


9.42 


1.6 


0.46 


0.4 


1.10 


3.47 


0.290 


NGC 129l( b > 


7.43 


9.45 


3.1 


0.77 


0.7 


0.70 


5.66 




NGC 1512 


7.25 


9.42 


3.2 


1.10 


1.0 


0.80 


0.87 


0.247 


NGC 1566 


8.21 


10.58 


4.6 


1.74 


1.5 


1.30 


1.13 


0.305 



a at 250 pm 

b with constraint f/ m j n > 0.7 
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TABLE 4 

Core flux density fractions for BLAST nearby 

GALAXIES 





Core 


Band 


Core 


Core Flux 


Galaxy 


Radius 


Flux 


Fraction 


['] [k P c] 


[urn] 


[Jy] 


[%] 



NGC 1097 


0.37 


1.8 


250 


13. 


35±0. 


03 


1 1 


93±0. 


.06 








oov 




Q7-t-0 


f)9 


12 




f)7 








ouu 


1 


ouznu 




o 




. uo 


NGG 1 991 


71 


2 


250 


I 


41i0 


02 


5 


42±0 


10 








350 





64±0. 


.01 


5. 


04±0. 


.17 








500 


0. 


23±0. 


.01 


3. 


58±0. 


.17 


NGC 1365 


0.29 


1.6 


250 


13. 


97±0. 


16 


9. 


58±0. 


.13 








350 


5. 


38±0. 


.10 


8. 


62±0. 


.19 








500 


1. 


35±0. 


.05 


5. 


.46±0. 


.23 


NGC 1512 


0.40 


1.2 


250 


1. 


,51±0. 


.01 


10. 


26±0. 


.08 








350 





,61±0. 


.01 


8. 


,79±0. 


.21 








500 


0. 


14±0. 


.01 


3. 


.43±0. 


.06 


NGC 1566 


0.98 


5.1 


250 


19. 


07±0. 


04 


37 


04±0. 


.19 








350 


7. 


.72±0. 


.04 


34. 


.78±0. 


.37 








500 


2 


,84±0. 


.02 


31. 


04±0. 


,16 


NGC 1808 


0.29 


0.9 


250 


10. 


86±0. 


01 


12 


94±0. 


.03 








350 


3. 


66±0. 


.01 


11 


06±0. 


.07 








500 


0. 


92±0. 


.01 


7. 


79±0. 


.05 



Note. — Flux uncertainties reported are measurement er- 
rors (o-(j) only. 



TABLE 5 

Resolved versus unresolved mass estimates for BLAST nearby 

GALAXIES 

Unresolved Mass' a ' Resolved Mass Ratio Aperture 

Galaxy . / j\/„ \ . / a/,. \ m,. Diameter 

log Ugj 1o h^J g n 

NGC 1097 8.37 9.17+°" 6.3j;^ 10.9 

NGC 1291 7.62 8 - 90 io.20 19 - 1+ -?'.o 13 - 2 

NGC 1512 7.25 8.38+.°,;i4 13-5±f? 5.8 

NGC 1566 8.21 8.92+°°l 5.1+];;? 6.8 



a DL07 model mass estimate from Table [3] 



